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Abstract We have succeeded in establishing a cosmo- 
logical model with a non-minimally coupled scalar field 
(j) that can account not only for the spatial periodicity 
or the picket-fence structure exhibited by the galaxy 
N-z relation of the 2dF survey but also for the spatial 
power spectrum of the cosmic microwave background 
radiation (CMB) temperature anisotropy observed by 
the WMAP satellite. The Hubble diagram of our model 
also compares well with the observation of Type la su- 
pernovae. The scalar field of our model universe starts 
from an extremely small value at around the nucle- 
osynthesis epoch, remains in that state for sufhciently 
long periods, allowing sufficient time for the CMB tem- 
perature anisotropy to form, and then starts to grow 
in magnitude at the redshift z of ~ 1, followed by 
a damping oscillation which is required to reproduce 
the observed picket-fence structure of the N-z relation. 
To realize such behavior of the scalar field, we have 
found it necessary to introduce a new form of poten- 
tial V{(j)) C!C (j)'^ exp{—q(f>'^), with q being a constant. 
Through this parameter q, we can control the epoch at 
which the scalar field starts growing. 
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1 Introduction 



It is currently believed, based on observations ( Riess et al.l 



1998HPerlmutter et al.lll999l : lRiess et al.l2004tlAstier et al 



i 



20061 iRiess et al.l 120071 ) of a number of Type la su- 
pernovae (hereafter referred to as SNIa), that the ex- 
pansion of our universe at present is in an acceler- 
ated stage. Furthermore, the measurements of the 
temperature anisotropy of the cosmic microwave back- 
ground radiation (hereafter CMB) by t he Wilkinson 



microwave anisotropy pro be (WMAP) (jSpergel et al. 



2007t iHinshaw et al.ll2007l ) have indicated that our uni- 
verse is flat. The high precision of these and other ob- 
servations have enabled us to intercompare and verify 
various candidate cosmological models. 

A flat ylCDM model consisting of Einstein's cosmo- 
logical constant A and cold dark matter, in addition 
to the ordinary baryonic matter and radiation is cur- 
rently the standard model because it can account rea- 
sonably well for the acceleration of the cosmic expan- 
sion as indicated by the SNIa's, the spatial power spec- 
trum of the CMB temperature anisotropy observed by 
WMAP satellite, and the primordial abundance of the 
light chemical elements produced through big-bang nu- 
cleosynthesis. 

Nevertheless, some key questions remain as to the va- 
lidity of the standard ylCDM model. For instance, the 
physical origin of the cosmological constant is presently 
unknown, although it might possibly be a component 
of vacuum energy. As such, we have yet to understand 
why the energy density of the cosmological constant re- 
mains unchanged with time. The numerical size of the 
yl-term estimated based on a particle physics consider- 
ation is generally m ore than 1 20 orders of magnitude 
too large (see, e.g., lRichll200ll p. 16) to be compared 



to the observed values. As an alternative approach 
to circumvent these problems of the A-teim, cosmo- 
logical models possessing quintessence or a phantom 



have been proposed (see, e.g.. iRatra fc Peebles 1988: 
Peebles fc R.atral Il988l : ICaldwelll 120021 : ICaldwell et al 



200 



3). These methods mterpret the dark energy in 



terms of the type of time- varying scalar field. 

The second consideration is that the cosmic age 
of the standard model, which is largely based on the 
WMAP 3-yr data, is supposed to be 13.7 Gyrs, whereas 
the possibility of it being, for instance, approximately 
14 ~ 15 Gyrs has been raised based on observations of 
certain globular clusters. We have already shown that 
we might be able to avoid this cosmic age problem by 
introducing a time-decaying constant as well as the ef- 
fect of g alaxy merger with a time scale of approximately 
3 Gyrs (JKomiva et al.ll2005l . i20061. However, in view of 
the various uncertainties involved in cosmic age deter- 
mination processes, 13 Gyrs appears to be a conserva- 
tive estimate for the lowe rbound for any cosm ological 
modeling, as indicated bv iKomiva et al.l (J2006I ). 

The third and most baffling problem, however, is 
the spatial periodicity of approximately 128/ij~qq Mpc 
with hioo = Ho km s-iMpc-V(100 km s^^Mpc^i) 
emerging in the relation of the galaxy number counts 
N and the redshifts z relation (hereafter referred to 
as the picket-fence s tructure of the N- z relat ion) orig- 
inally discovered by iBroadhurst et al.l (jl990l ) through 
a pencil-beam survey of field galaxies. In each of the 
four directions of the survey, an epoch of relatively high 
galaxy number counts and another that exhibits com- 
paratively low number counts appear alternately, as 
though concentric spherical shells (or epochs) of higher 
galaxy number densities were surrounding us with their 
individual centers situated at our location. 

Such an observed characteristic feature, if it were re- 
flecting an actual physical concentration of galaxies ex- 
isting at certain distances from us, would deflnitely be 
incompatible with the cosmological principle that pre- 
sumes the uniformity and isotropy of our space-time. 
In fact, based on a large number of numerical simula- 
tions using the Einstein-de Sitt e r mod el as well as the 
ylCDM models, lYoshida et al.l (|2001[ ). among others, 
demonstrated that the probability to obtain a periodic 
structure in the spatial distribution of gala xies such 



as that inferred from the results of Broadhurst et al. 



( 19901 ) should be less than 1/1000. Hence, it is often 
the case tha t the existence of t he picket-fence structure 
reported by IBroadhurst et al.l (|1990l) is dismissed as an 
observational fluke. 

However, this problem is not so simple. The N-z re- 
lation obtained from the 2dF Galaxy Redshift Survey 
(2dF GRS) also exhibits a definite picket-fence struc- 
ture at th e interval of rough ly Az = 0.03 (see, e.g., 
Fig. 17 of lColless eraDboOlh . The locations of higher 
galaxy number counts are found at 2; ~ 0.03, 0.06, 0.08, 



0.1 1, . . . , which is essentia lly the same as those found 
bv IBroadhurst et al.l (1990). Even in the N-z relation 
obtained through the Sloan Digital Sky Survey (here- 
after SDSS), protrusions still appear in the locations of 
relatively higher galaxy number counts compared with 
th ose predicted from th e standard model (see Fig. 2 
of iTegmark et al.l (|2004l )) at z ~ 0.03, 0.05, 0.08, 0.11, 
. . . , although their amplitudes of variation are much 
less pronounced compared to those observed in the 2dF 
GRS data. 

For the quantitative comparison of our theoretical 
computations, we would like to adopt the N-z rela- 
tion taken from the 2dF GRS data for the following 
reasons: i) the sky coverage of the 2dF data is more 
symmetric with respect to the celestial equator than 
that of the SDSS data, including not only a major 
portion of the Southern hemisphere sky, but also a 
comparable portion of the Northern hemisphere sky, 
w hereas the renderings of the SPS S project described 
bv lAdelman-McCarthv et al.l (|2007l l indicates a notable 
asymmetry in the hemispheric coverages, and ii)the 
width of the bin for number count statistics used for 
the 2dF data is narrower by a factor of approximately 
2 compared with that for the SDSS data. 

One viable explanation for the picket-fence struc- 
ture of the N-z relation within the framework of 
the Robertson- Walker metric would be the intro- 
duction of a scalar field that non-minimal ly couples 
with t he cu rvature scalar, as proposed by iMorikawal 
(|1990l I199I ). Morikawa's model was su bsequently ex- 
tended by Kashino fc Kawabatal (|l9941) . and later by 



Fukuvama et al.l (|l997l) . to include the effects of mat- 
ter and the A term. Such coupling can produce the 
picket-fence s tructure in the followin g way: the scalar 
field model of lMorikawal ( 1990Lll99ll ) assumes a poten- 
tial V{(j)) proportional to (j)"^ so that the state of our 
universe oscillates back and forth around the potential 
minimum, inducing an oscillation of the Hubble param- 
eter as well, due to the curvature coupling of 0. Hence, 
an epoch of relatively rapid expansion and an epoch of 
relatively slower expansion appear alternately due to 
this oscillation. During the former epoch, the number 
density of galaxies would diminish noticeably; whereas 
during the latter epoch, the decrease in the number den- 
sity would remain less prominent. We would then get 
a picket-fence structure as an apparent or illusionary 
effect in the N-z relation, because we are observing the 
two types of epochs alternately with increasing value of 
z. 

Since the spatial distribution of galaxies remains 
homogeneous for a given instance of time, the cos- 
mological principle can still be preserved. Note that 
iFukuyama et al.. (jl997ii) pointed out that the adoption 
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of negative values for the coupling constant ^ yields a 
succession of mini-inflationary (concave shaped) cosmic 
expansion with time, thereby prolonging the cosmolog- 
ical age. We shall therefore employ negative values for 
^. The temporal oscillations of the fundamental con- 
stants, such as the gravitational constant G, th at would 
be caused by the scalar field (,Hill ct al.lll99Cir ) and the 
possible presence of a coherent peculiar vel ocity field 



are l ikely to enhance that apparent effect (JHill et al. 
199ll ). 



On the other hand, based on the pri mordial abun 
dance s of the li ght chemical e l ement s, ISalgado et al 



(|l996l ). and later lOuevedo et alj (|l997[ ). argued that the 
scalar field needs to be at an extremely low level dur- 
ing the time periods, including that for the big bang 
nucleosynthesis, and demonstrated that under certain 
conditions, it is possible to suppress the scalar field (j) 
to an extremely small value {<j) < 10^^") over certain 
periods of time at the approximate time of nucleosyn- 
thesis, which subsequently starts to grow in magnitude 
and simultaneously exhibits damping oscillations as we 
move toward the present. 

We shall therefore restrict ourselves to the models in 
which the scalar fields undergo stationary states with 
sufficiently small values over the relevant periods of 
time. One notable problem associated with comput- 
ing such cosmological models starting from the present 
epoch, as is usually done, is that one must perforin a 
fine-tuning for the initial value for cj) io a, large number 
of decimal places in order to arrive at the stationary 
state of (j) at or near the time of big bang nucleosynthe- 
sis. However, we have found that by solving the evo- 
lution equations from the time of this nucleosynthesis, 
when (/) ~ 0, toward the present, we need not carry out 



any su ch fine-tuning, as was suggested bv lSalgado et al 
( 19961 ). The theoretical basis and a number of sample 



computations will be published elsewhere. 

In general, we cannot hope to sufficiently constrain 
the parameter values unless more than one type of ob- 
servational data are employed, as we ourselves have 



learned in other investigations (jKomiva et al.l l2006f ). 
Therefore, we also make use of another crucial obser- 
vation, viz., the spatial power spectrum of the CMB 
temp erature anisotropy obtained from the WM AP 3-yr 
data (jSpergel et al.ll2007t iHinshaw et al.ll2007n to nar- 
row the possible ranges of the model parameters. The 
candidate model thus obtained will be further cross- 
checked against the observed Hubbl e diagram of SNIa 's 
produced with the Gold data set of iRiess et al.l (|2007l ). 
The purpose of the present study is thus twofold. 
First, we investigate whether it is possible to determine 
a cosmological model possessing the non-minimally cou- 
pled Morikawa scalar field that can account for both the 



picket-fence structure of the N-z relation of the field 
galaxies and the power spectrum of the CMB tempera- 
ture anisotropy given by the WMAP 3-yr data. Second, 
we attempt to determine a simple method by which to 
modify the basic equation or equations to construct a 
model that can better serve our purpose, if the above- 
mentioned model is found to be inappropriate. We 
herein demonstrate that the use of a potential for the 
scalar fiel d that is somew hat different from that em- 
ployed bv iMorikawal (|l99d ) is better for this purpose. 

The remainder of the present paper is organized as 
follows. In Section 2, we present the set of evolution 
equations for our universe, incorporating the newly in- 
troduced potential for cf), and discuss the characteristic 
features of our cosmological model compared to those 
of the standard /ICDM model. In Section 3, we in- 
vestigate the effect of each of the principal model pa- 
rameters, including the power index q employed for our 
newly proposed potential on the theoretical computa- 
tions of the N-z relation and the power spectrum of 
the CMB temperature anisotropy, and show that there 
exists a set of parameter values with which we can re- 
produce, rather satisfactorily, both the 2dF data and 
the WMAP 3-yr data. In addition, we demonstrate 
that the model thus obtained can account for the ob- 
served Hubble diagram of the recently observed SNIa's, 
at least to a degree comparable to that of the standard 
model. In Section 4, we summarize our results and 
present our conclusions. 



2 Basic Equations for Theoretical 
Computations 

2.1 Derivations of Time Evolution Equations 

The relevant Lagrangian L for our cosmological model 
is as shown below: 



L ^ 



1 /'m^,c\'^ 



R 



1 



g'^-'d^i^d.^p 



^(mr) 



47rG 



■v-^ 



SttG 



zA 



(1) 



where g^'^ are the (/i, :/)-components of the metric ten- 
sor, ij} is the scalar field, R is the scalar curvature, <^ 
is the coupling constant, 5 is a scalar constant, m^ is 
the mass of the scalar field, G is Newton's gravitational 
constant, c is the velocity of light, A is Einstein's cos- 
mological constant, and L^™"^) is the Lagrangian due 
to matter and radiation. In addition, we specify that 
d^-ip = dip/dx^. Each Greek subscript ranges from 



to 3, and Einstein's summation convention is implicitly 
assumed in the above equation, as well as the equa- 
tions presented hereafter. Then, based on the principle 
of least action, we extremize the action integral 



S- 



d x\/—g L 



(2) 



by setting 6S = 0, which leads us to the following grav- 
itational field equations: 



&t:G 



d^^ipd^ip - -^giiudptpOPip 






T^(mr) 



(3) 



where the symbol "D" signifies the d'Alembertian, the 
symbol "(• • • );;i" is the covariant derivative of a quan- 
tity "•••" with respect to spatial coordinate a;^, and 
the quantities T^™"^ are the (jjl, i^)-components of the 
energy-momentum tensor of matter and radiation. 

Note that the potential term (the fourth term on the 
right-hand side of Eq.(IT])) refl ects our niqdification t o 
the original form employed by iMorikawal ( 199CI Il991l ). 
which coincides with the latter if q = 0. As we shall 
show later herein, this new form of potential enables 
us to control the epoch when the growth of the scalar 
fleld -0 begins to take place. The closer this epoch is 
to the present era, the less affected are the amplitudes 
of the spatial power spectrum of the CMB tempera- 
ture anisotropy in the large-scale domain because of 
the reduced effect of t he late-time integrat ed Sachs- 



Wolfe effect 
on 



see. 



.Sii 



NaselskvetailboOGl ). Based 



Kashino fc Kawabatal (|1994| ). we have included the 



A term, viz., the flfth term on the right-hand side of 
Eq.Ill) because it provides greater flexibility, yields a 
better fit to the CMB spectrum, and makes the age 
of our universe older by dominating the time periods 
prior to the resurgence of the scalar field (/). The fact 
that we have a stationary state of <f> at an extremely 
small value does not necessarily mean that the universe 
was born with ~ 0. On the contrary, the Big Bang 
is supposed to have been accompanied by a scalar field 
of considerable magnit ude, as has been pointed out by 
Fukuvama et all (|l997h . 

In accordance with the inflationary cosmology, we 
restrict ourselves to a flat geometry, so that the 



Robertson- Walker line element ds takes the following 
form: 

ds^ = c^df ~ a^{t)[dx^ + x^ldei^ + sin^ 9 d(t>'^)]. (4) 

Here, a(t) is the cosmic scale factor at proper time t 
normalized to the present-day value, whereas Xj ^j and 
(j) are the comoving angular coordinates. The matter 
included in our model universe is assumed to consist 
of baryonic matter as well as cold dark matter (CDM), 
but not massive neutrinos. Furthermore, the scalar field 
is assumed to interact with the gravitation merely by 
coupling with the scalar curvature, so that there is no 
mutual energy exchange between the matter, radiation, 
and scalar field. Under this assumption, the equation 
of state for each constituent is given as follows: 



p;„c^ = 3-(/9„c^ -t- P,„) = 0, 
a 

PrC — 3—{prC +Pr)—0, 



P^C 



3-(p^c2 + P^) = 0, 



together with 



rjYi — U, jV — t^PrC 



(5) 
(6) 

(7) 
(8) 



where the quantities p and P stand for the mass den- 
sity and the pressure of each constituent of the uni- 
verse, respectively, and the subscripts m, r, and ip 
designate "matter" (baryon-|-CDM), "radiation", and 
"scalar field", respectively. We also employ a short- 
hand notation A = dA/dtc, where ic is the conformal 
time defined as 



, c , , ■ dA adA 

tr — I —dt, so that A — - — = — — . 
a dtr c dt 



(9) 



On the other hand, the mass density p^ and the 
pressure P^ for the scalar field tp resulting from the 
principle of least action are as shown below: 



p,pc 



2 a^ a* a-^ 



-im'*' 



exp 



AnG 



ri (10) 



Pi, 



2 a^ a 



3 



a^ 



-2^J_^^_2^^(V-)2-2e:^VV' 



a^ 



O'^ 



-\m'*''A-^'-^*'} 



(11) 
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Substituting Eqs.(ITO|) and ([TT]) into Eq.®, we obtain 
the time evolution equation for the scalar field: 

a a \ n / 

(12) 

where we have employed a dimensionless quantity 4>, 
instead oi ip: 



AtiG 
"3^ 



4,. 



(13) 



From the (O,0)-components and (l,l)-components of 
Eq.Q, and using Eg. p^ . we obtain the following set of 
the time evolution equations for the cosmic scale factor 
a: 



6C# + {(6^00)2 + (1- 6^0') 



— 5" 1 5- + "/l.oa 



1/2 



/(I -6^02), 



(14) 



a 
a 



.(1 _ 6^0") - 3(<^)2 - 2A^-^^ + 6^(0)' 



-eCa^ (^)'(l -#')</>' exp (-#2) 

~2~^ \ a ~^ 

/{I -6^0^(1 -60}, (15) 

where Hq is the Hubble constant. The quantities ^ni,o, 
flr.Oi ^-nd rjyi.o are the present mass densities of mat- 
ter, radiation, and the cosmological constant A nor- 
malized, respectively, to the present value of the crit- 
ical density Pc,o(= iH^/SnG), viz., flm^o = Pm,o/pc,o, 
ftr,o = Pt,o/Pc,o and flA,o — {c/Hq)'^A/3. The density 
parameter of the scalar field "0 for the present epoch is 
similarly defined as Sl^^o = Pv.o/Pc.o- Since the mass 
density p^ in terms of cj) should be equal to p^ , we have 
^4>,o = ^i/),o- Furthermore, since we are concerned only 
with the flat geometry case, we obtain the following 
constraint from Eq. (|14p for the present era: 



^m,0 + ^r,0 + ^0,0 + ^yl,0 — l, 

SO that 



(16) 



(17) 



which is to be substituted to Eg. ([14] 



2.2 Stationary State and Onset of Damping 
Oscillation of Scalar Field 



Salgado et all (|l996l ) and lOuevedo etldl (|l997f ) indi- 
cated that the magnitude of the scalar field must re- 
main extremely small during the epochs around the big 
bang nucleosynthesis in order not to spoil the successful 
prediction made on the cosmic abundances of light ele- 
ments based on the Einstein-de Sitter model universe. 
However, to give rise to such a plateau state of the 
scalar field (j) at those periods of time without under- 
taking any extreme fine-tuning of the presentday value 
of (j) n such as those carried out in Salga do et al.l (jl996l ) 



and lOuevedo et al.l (|1997I ). we would like to start inte- 
grating the evolution equations at the time ic,init corre- 
sponding to a = 10~^, which is sufficiently close to the 
big bang nucleosynthesis toward the present time ic,o- 

In view of the fact that the universe at this stage is 
dominated by radiation, we approximate the values of 
a(ic,init) and d(ic,init) with those given by the ylCDM 
model for a given set of the values of flm.o, ^r.o, and 
flA,o- Furthermore, we take 0(ic,init) — because of 
the assumed stationary state of (j). Then, (/)(tc,init) is 
the only parameter that needs to be specified at tc.init, 
which is determined in the following manner: a starting 
value for 0(tc,init) is picked, and the model computation 
is carried out all the way to the present time, and the 
resulting value of the Hubble parameter Hq is inspected 
to see if it agrees with 72 km s~"'^Mpc~^, the intended 
value, within 0.1%. If not, we multiply (/)(ic,init) by 10, 
and perform the model computation again. This pro- 
cedure is repeated until we bracket 72 km s~-^Mpc~^, 
after which the Newton method approximating Hq — 72 
with a piece- wise linear function of log 4>init is applied 
to find the desired value of (?!)(ic,init)- 

A sample behavior of the scalar field (f> is shown 
in Fig. 1 as a function of scale factor a for a model 
universe with Morikawa scalar field inodc\{q = 0) ob- 
tained using the procedure described above: we have 
employed f^,n,o = 0.237, f^^.o = 8.066 x 10"^ corre- 
sponding to TcMB,o = 2.726K, ^ — —40, and mg — 
3.7 X 10"'^-^ft-ioo eV. Note that cosmic age ijnit — 23.86 
s (yielding ic,init = 4.63 x 10"^ Mpc), d(ic,init) = 
2.157xl0-^and(/)(tc,init) = 2.341 x 10"^" at a = 10~^ 
from which Hq = 72.031 km s^^Mpc^^ and 0o = 
7.979 X 10"^ resuh. 

Next, we want to demonstrate that the time for (f> to 
start growing can be controlled by changing the value 
of the power index q introduced in Eq. (fn| . Figure [U^a) 
shows the variation of (P as a function of scale factors 




Fig. 1 An example of the variation of tiie Morikawa-type 
scalar field as a function of scale factor a, starting from 



a stationary state at a 
description) . 



10 (see the text for the model 



a for three values of q, viz., 0, 6, 000, and 12, 000. The 
values of the model parameters employed are as fol- 



n„ 



0.237, n^^o = 0.15(hence fl 



A,0 



0.613), 



lows: 

hioo = 0.72, ^ = -40, and m^ = 3.7 x lO-^'^/iioo eV. 
They are taken from our best-candidate model deduced 
based on a large number of trial and error search in com- 
parison with the observational data of the galaxy N-z 
relation and the spatial power spectrum of the CMB, 
which we shall discuss shortly. Note that the onset of 
the damping oscillation of 0^ is delayed as the value of 
q increases. Figsj^Ib), [H^c), and^fd) indicate the cor- 
responding derivative of with respect to the Hubble 
time Hot, Hubble parameter H/Ho{see also Figl^l an 
close-up view of this diagram), and effective gravita- 
tional constant GcS given by 



Gc« 
G 



1 



1 - 6^02 



(18) 



respec tively, where G is Newton's gravitational con- 
stant ( Morikawalll991 1^l. These quantities exhibit sim- 
ilar oscillatory behaviors as functions of a to that of 

W,2 



^ The scalar fi eld ip used in the present work corresponds to <j> 
in iMorikawal l|l990l . ll99U ). Hence, the quantity used in the 
present study is equal to \/47rG/3c^ times Morikawa's </». See 
also Eq. H13|l of the present work. 
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Fig. 2 (a) Dependence of 0^ on the value of q. The 
abscissa is the scale factor a; (b) Dependence on q of the 
derivative of (p with respect to Hubble time Hot. The ab- 
scissa is the scale factor a; (c) Dependence on q of the onset 
time of the oscillation of the Hubble parameter H/Hq nor- 
malized to the Hubble constant. The abscissa is the scale 
factor a; (d) Dependence on q of the time for the onset 
of oscillation of the effective gravitational constant Gcs/G 
normalized to the Newton's gravitational constant G. The 
abscissa is the scale factor a. 
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4; 




Fig. 3 Detail of the behavior of H/Hq depicted in Figl^Jc) 
near the recent epochs. 



3 Observational Tests 

3.1 Galaxy Number Counts with Respect to Redshift 

As has been mentioned ear her, the observed data on 
the galaxy number counts as a function of redshift z 
obtained through t he 2dF survey (see. e.g., Fig. 17 of 
Colless et al.l 120011 ) exhibit alternating regions of com- 
paratively high counts and comparatively low counts at 
the interval of approximately Az = 0. 03, essentially the 



same feature as that discovered by iBroadhurst et al 



(1990|) through their pencil-beam survey, which is re- 
ferred to herein as the picket-fence structure. It should 
be emphasized that the standard model, namely the flat 
ylCDM model derived on the basis of the WMAP obser- 
vations, is not compatible, or is at least extremely diffi- 
cult to reconcile with this picket-fence structure, unless 
we are willing to set aside the cosmological principle. 

In this section, we compare the N-z relationships for 
the bj-band computed using the proposed model having 
the new p otential for th e scalar field and the 2dF sur- 
vey data ( Colless et al.ll2G01 ). The theoretical compu- 
tations are carr ied out using the computer code ncmod 
( Gardneiill998r ) with a modification to incorporate our 
cosmological model. Briefly speaking , we adopt the fol- 
lowing valu es for the Schechtei J 19761 ) parameterization 
derived by iNorberg et al.l ([2002) based on the fitting 
to the 2dF data: M^^ - 51ogio/iioo = -19.66, a = 
-1.21, $* = 1.61 X W'^hlooMpc'^. Here, it is as- 
sumed that galaxies fainter than rrib, = 19.45 cannot be 
seen, which is taken into account in our computations 
in terms of the selection function. Furthermore, the 
spectral energy distributions (SEP ) of different types 
of gal axies are those employed by lYoshii fc Takahara 
( 19881 ). 
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Fig. 4 Dependence on the coupling constant ^ of the com- 
puted N-z relation of galaxies. The abscissa indicates the 
redshift z, while the ordinate is the fractional number count 
A''. For comparison, the 2dF data are also plotted. 

First, let us investigate the dependence of the theo- 
retical N-z relationship on each of the model param- 
eters. Figure [4] shows the theoretical N-z relation- 
ships computed for ^ = 0, —20, —40, and —80 us- 
ing the proposed model with the new potential. The 
values of the other model parameters are as follows: 
fim^o = 0.237, n^fi = 0.15 (hence ^a,o = 0.613), 
hioo = 0.72, TO^ = 3.7 X IQ-^'^hwo eV, and q = 12, 000, 
which are taken from our best-fit model (see FigfTO|. 
which is capable of reproducing both the 2dF data and 
the WMAP 3-yr data. Note that no spatial periodicity 
arises for ^ = 0, but that the amplitude of variation 
of the galaxy number count N with z tends to be en- 
hanced as the absolute value of ^ is increased. In ad- 
dition, a slight change in the period of oscillation takes 
place with the value of ^ as well. As a result, ^ = — 40 
is selected. 

Figure m shows the dependence of the N-z relation- 
ship on the density parameter of the scalar field ^^fi, 
for which we have employed 0.00, 0.15, 0.30, and 0.60 
together with ^ = —40, whereas the values of the re- 
maining parameters are kept the same as those em- 
ployed for Fig. 21 Note that the amplitude of variation 
of N increases with increasing fi^^o- The case in which 
^4>,o = is found to be essentially the same as that for 
the flat ylCDM model with fim.o = 0.237 because the 
value of remains sufficiently small throughout the rel- 
evant epochs. We select V,^^ = 0.15 as the best choice. 

Figure [6] shows the theoretical N-z relationships for 
TO0 = 2.0 X lO-31/iioo eV, 3.7 x lO-^i/iioo eV, and 5.4 x 
10~^-^/iioo eV, respectively. The values of the remaining 
parameters are the same as those for Fig. [31 except that 
we have applied l^^.o = 0.15. The period of oscillation 



6 


1 l| 


— r 


1 
1 

1 


1 1 1 

2dF GRS 


b 

I 


.'a 
n'l 


\ 


fl 


n^.o = 0.00 

%,0 = 0.15 - 

%o = 0.30 


X 


M 


fi 


% 


"#.0 = 0.60 


1 

o 

X 


-/ 1 


1 


7 

o 


/ 1 1 1 


\_ " 



2dF GRS 
q» = 
g = 6,000 
g = 13,000 



0-05 



0.1 



0,15 



0.2 



0.25 



0.3 



0.35 




0.35 



Fig. 5 Dependence on £7^,0 of the computed N-z rela- 
tionship of galaxies. The abscissa indicates the redshift z, 
while the ordinate is the fractional number count A^. For 
comparison, the 2dF data arc also plotted. 



Fig. 7 Dependence on the power index q of the computed 
N-z relationship of galaxies. The abscissa indicates the red- 
shift z, while the ordinate is the fractional number count A'^. 
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Fig. 6 Dependence on the mass m^ of the scalar field, 
of the computed N-z relationship of galaxies. The abscissa 
indicates the redshift z, while the ordinate is the fractional 
number count A''. For comparison, the 2dF data are also 
plotted. 



of N with respect to the redshift z tends to decrease 
and its amplitude becomes smaller with increasing m^. 
We select 3.7 x 10~'^i/iioo as the best value for m^. 

Figure [7] depicts the theoretical N-z relationships 
computed for q = 0, 6,000, and 12,000, respectively. 
The values of the other parameters are the same as 
those used for FigEl except that m^ = 3.7 x 10~^^hioo 
is used. Note that the amplitude of variation of N is 
only slightly affected by the value of q, although the pe- 
riod of oscillation changes slightly. In other words, we 
could maintain the amplitudes of N variation at suffi- 



cient magnitudes even if we vary the value of q over an 
appreciable range. Nevertheless, in the next section, 
we shall show that the value of q significantly affects 
the large-scale portion of the spatial power spectrum of 
the CMB temperature anisotropy, which enables us to 
select 12,000 as the best value for this parameter, and 
arrive at a model that can account for both the picket- 
fence structure of the N-z relationship of galaxies and 
the observed power spectrum of the CMB temperature 
anisotropy. 

Note that the ylCDM model predicts only a smooth 
N-z curve without abrupt changes, just like the 
case of ilcjifi = in FiglS] Unlike the cases for 
the original M orikawa models ( Morikawal Il990l . Il991 ; 
Fukuyama et a l. 1997), the scalar field in the present 



models begins to exhibit significant oscillation of appre- 
ciable amplitude at approximately z ~ 1, but is rapidly 
damped toward the present. Note that as the value of 
q increases, the onset of the oscillation shifts toward 
a more recent epoch. However, since we observe the 
galaxies only up to z = 0.35 or so in the 2dF data, 
an exact choice of the q value is not important with 
respect to the theoretical N-z relationship. However, 
later herein, we shall demonstrate the fundamental im- 
portance of the non-zero q in simultaneously producing 
any satisfactory fit of the theoretical spectrum to the 
WMAP data on the CMB temperature anisotropy. 

3.2 Cosmic Microwave Background Radiation 

The high-resolution measurement of the spatial power 
spectrum of the temperature anisotropy of the CMB by 
the Wilkinson microwave anisotropy probe (WMAP) 
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should provide us with another crucial test for the pro 
posed model. For t his purpose, we emp l oy the WMAP 



npi 

m 



Hinshaw et al. 



3-yr project data (jSpergel et al.1 1200 
2007t ). 

For the set of our cosmological models that have 
provided reasonable fits to the observed galaxy N-z 
relation, we calculate the theoretical angular power 
spectrum of the CMB temperature anisotropy, while 
making use of a modified Boltzmann code CAMB 
( Lewis et alJ 120001 ) . which is based on a line-of-sight 
integ ration approach, as employ ed in the CMBFAST 
code (jSeliak fc Zaldarriagalll996l ). We first tabulate the 
values of {a, a, a, cj), 4>, and 0} obtained through in- 
tegrating Eqs. p^ . p^ . and pS]) with respect to the 
conformal time. These tables are then used to perform 
interpolations to find the values of these variables for 
a given time inside the power spectrum computation 
program. 

As for the perturbation compu tations, we make 
use of t he formulation generated bv iBaccigalupi et al 



|2000l ): IPerrotta fc Baccigalupil ([200^. The required 
perturbation equations can be obtained by substitut 



mg 



F = 



SttG 



^r 



_ 1 /■ m^C 
~ 2 



tp exp 



47rG 2 



(19) 



(20) 



in their formulae (see Appendix). We have also rewrit- 
ten the part of the perturbation equations in the CAMB 
usi ng the quantities appr o priate for our model. As 



Perrotta fc Baccigalupil (|l999f ). the adiabatic con- 



dition is thereby assumed as the initial condition for 
perturbations. Since the scalar field is practically non- 
existent in the early epochs of the proposed models, 
the effect of (f> can be ignored in comparison with the 
other components comprising the proposed models, so 
that the universe can be represented with nothing more 
than a radiation-dominated ylCDM model, as is shown 
in FiglH which indicates the fractional contribution of 
matter (baryonic matter -I- CDM), radiation, yl-term, 
and scalar field (j) to the total density parameter, which 
is unity because we are herein concerned only with the 
flat geometry. 

This diagram shows that radiation is dominant over 
the other constituents for a < 3.4 x lO""*. The val- 
ues of the model parameters employed here are as fol- 
lows: i^hfihloo = 0-024, Qn^fl = 0.237, Hq = 72 km 
s-^Mpc-'i, r^^^o = 0.15, nA,o = 0.613, ^ = -40, 
m^ = 3.7 X 10~^^/iioo eV, and q = 12,000, which char- 
acterize the best-fit model proposed herein. Hence, the 
initial perturbations of the densities of the baryonic 



matter, CDM, and radiation, as well as those of the 
metrics, m ay well be assumed t o be given with those 
derived by iMa fc Bertschingeil (1995") for the ylCDM 
model. In what follows, we shall discuss the theoreti- 
cal power spectra of the CMB temperature anisotropy 
computed for the proposed cosmological model. 




Fig. 8 Fractional contribution of radiation, matter, A- 
terni, and scalar field (f) to the total density parameter at 
each epoch. The total density parameter is always unity 
because of the flat geometry adopted herein. 



3.2.1 Cosmological model with Morikawa's scalar 
field: q = case 

Next, let us consider whether we can select an ade- 
quate model or models with Morikawa's scalar field 
from among those that can reproduce the 2dF GRS 
data based on the comp arison with the WMAP 3-yr 
data dHinshaw et al.l 120071 ). The results of experiments 



along this line ( Kawabata et al.l l2002t iHirano et al. 



I2OO6) indicated that it is extremely difficult, if not im- 
possible, to come across such a model because the pres- 
ence of a scalar field that is large enough to account for 
the picket-fence structure tends to yield too large a tem- 
perature anisotropy of the CMB in the large-scale do- 
main, as shown in Figl9l Here, we employ ^ = — 1 and 
—3, the absolute values of which are an order of magni- 
tude smaller than that inferred from the previous analy- 
sis of the N-z relation of galaxies (^ = —40). The values 
of the other parameters are as follows: fini,o = 0.237, 
^hfihloo = 0.024, 17^,0 = 0.15(hence fl^.n = 0.613) , 



m^ 



3.7x10 



-31; 



eV, h 



100 



0.72 ( Freedman et al 



20011 ). ris = 1 is the spectral index (Harrison-Zel'dovich 



spectrum) for the initial condition for th e matter den- 
sity perturbation, Tcmb.o = 2.726 K (JMather et al 



10 

Il999) for the present CMB temperature, and the op- 
tical depth of our u niverse is r = 0.089 for the CMB, 
(jSpergel et alj 120071 ). which, together with ^ = —40, 
characterize the best-fit model proposed in the present 
study. Furthermore, we adopt Yhc = 0.24 for the cos- 
mic abundance (by mass) of He and 3.04 for the number 
of neutrino species (JNagata et al.ll2004l ). 
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Fig. 9 Power spectra of the CMB temperature anisotropy 
computed for the Morikawa scalar field model with ^ — —1.0 
and —3.0. For the purpose of comparison, the theoretical re- 
sult obtained by Spcrgol ct al. ( 200*^) for their ^ CDM model 
and the WMAP 3-yr data lJHinshaw et ahllJOOTh on which it 
is based are also shown. 

Also shown in FiglH] is the theoretical spectrum 
computed for the ylCDM model adopting r^b.o^ioo = 



0.02229, n 
0.958, and r 



oh^ 



100 



0.1277, h 



100 



0.732, 



0.089 derived by the WMAP team 



( Spergel et all l2007^ together with the WMAP 3-yr 



data. Note that the locations of the peaks of the spec- 
trum computed for the model with the Morikawa scalar 
field tend to be shifted in the direction of the larger 
scale side in comparison with those computed for the 
ylCDM model. The reason for this is that the horizon 
scale of the scalar field model at the time of recom- 
bination is greater compared with that of the ylCDM 
model because the effective gravitational constant Goff 
is smaller in the past. T he opposite is t rue in the case of 
the scalar field model of lNagata et al.l ( 20041 ) . for which 
the value of GeS is larger in the past. 

Note also that the power spectra obtained for this 
scalar field model are significantly larger in the large- 
scale domain than that obtained for the ylCDM model. 
This may be because the equation of state for the scalar 
field models fluctuates with time, so that the late-time 
Sachs- Wolfe effect is enhanced. 



Hence, we are inevitably led to the conclusion that 
we cannot possibly hope to reconcile the existence 
of the picket-fence structure in the N-z relation and 
the spatial power spectrum of the CMB temperature 
anisotropy within the framework of the cosmological 
model possessing the Morikawa scalar field. 

3.2.2 Cosmological model with a new potential for the 
scalar field: q =/= case 

Note that we have made the scalar field stationary at 
nearly null level for a roughly less than 5 x 10~^, even 
in the case of g = (see FigllJ. It may therefore be in- 
ferred that the magnitude of the scalar field must still 
be too large during the formation of the spatial power 
spectrum of the CMB. In other words, we might be able 
to obtain a closer fit to the WMAP data in the large- 
scale region, if we could retard the growth of the scalar 
field (j). Based on theoretical considerations and nu- 
merical verifications, we developed a simple means by 
which to achieve this objective, namely, the introduc- 
tion of a new form of potential V{(f>) for the scalar field, 
which is proportional to (j)'^ exp (— q(^^) with q being con - 
stant instead of 0^ adopted bv iMorikawal (|l990lll99ll) . 
as has been indicated in Eqs.pl]) and UH). Moreover, 
our potential tends to Morikawa's as q approaches 0. 

Recall that the different choices for the values of q 
hardly affect the resulting theoretical curves of the N- 
z relation, as we have demonstrated in Fig. [T] The 
actual process we used to obtain the final candidate 
model is a sort of trial-and-error approach in which the 
quality of fit of the resulting N-z curve and the power 
spectrum of the CMB to the observational data are si- 
multaneously inspected. Rather than delineating the 
technical details here, however, we simply adopt the 
best-fit model that has been obtained and carry out 
a sensitivity study for each of the model parameters 
(g, ^, ri^^Oi JTi^j ^m,Oj and f^b.o) in order to illustrate 
how we can select a particular value for each parame- 
ter. Therefore, our procedure does not guarantee that 
the true best-fit model, if any, can be located. This 
does not matter, however, because the main objective 
of the present work is to demonstrate that a candidate 
model or models that account for both the picket-fence 
structure and the observed temperature anisotropy of 
the CMB can be found in a unified manner. 

We shall employ our best-fit model as the basis for 
our subsequent analysis, as was done in the 9 = case. 

First, let us vary the value of q and compute theoret- 
ical power spectra of the CMB temperature anisotropy. 
The results are shown in Fig. [TUlfor 9 = 6,000, 9,000, 
and 12,000. Also shown are the be st-fit curve from 
the ylCDM model |Spergel et al.ll2007t l and the WMAP 
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3-yr data (H inshaw ct a l. 2007) . We are clearly able 
to produce CMB spatial power spectra, even with the 
scalar field model, that ar e at le ast as good as those 
obtained by ISpergel et alj (|2007( ). Hence, we adopt 
q = 12,000. However, the parameters q, flh,Oi f^m,0: 
and the optical depth r only slightly influence the am- 
plitudes and spatial periods of the picket-fence struc- 
ture of the N-z relation. 
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Fig. 10 Spatial power spectra of the CMB tempera- 
ture anisotropy computed for our scalar field models with 
q =6,000, 9,000, and 12,000, respective ly. For the purpose o f 
comparison, the spectrum obtained bv lSpergel et al.l ([20071) 
for the ACDM model as well as the WMAP 3-yr data are 
also reproduced. 



Let us next investigate the sensitivity of the coupling 
constant ^ to the shape of the spatial power spectrum 
of the CMB temperature anisotropy. 

FigurcfTTJ shows the theoretical spectra computed for 
£, ^ -20, -40, -60, and -80, but with q = 12, 000 for 
all of them. The values of the other parameters are 
the same as those employed for Fig. [TOl The WMAP 
3-yr data are also shown in the diagram for compari- 
son. The theoretical spectra with ^ = —40 and —20 are 
now in good agreement with the WMAP data. This is 
in contrast to the case with q = 0, where ^ = —3, or 
even —1, gives rise to an appreciable deviation from the 
observational data in the large-scale portion, presum- 
ably due to the late-time Sachs- Wolfe effect. Taking 
into consideration the resulting N-z relation, we adopt 
^ — —40 rather than —20, as indicated in Fig. ID 

Figure [T^ shows the theoretical power spectra of the 
CMB temperature anisotropy for Q^^o — 0, 0.15, 0.30, 
and 0.45, respectively. The values of the remaining pa- 
rameters are the same as those used for Fig. [TTl except 
that we have employed ^ = —40 here. Note that the 
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Fig. 11 Spatial power spectra of the CMB temperature 
anisotropy obtained for ^ = —20, —40, —60, and —80. The 
values of the remaining parameters are the same as those 
employed for Fig. 1101 except that q — 12, 000 is used here. 
Also shown for comparison are the WMAP 3-yr data. 
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Fig. 12 Dependence on 57,^,0 of the spatial power spec- 
trum of the CMB temperature anisotropy computed with 
our model. The values of the other parameters are the same 
as those employed for Fig. 1111 except that ^ = —40 is used 
here. For the purpose of comparison, the WMAP 3-yr data 
are also shown in the diagram. 



large-scale portion of the computed spectrum depends 
rather critically on the value of the density parameter 



fltj)fl of the scalar field. We thus select Q^^ = 0.15 as 
the most appropriate value. 

Figure [13] shows the CMB power spectra obtained 
for m^ = 2.0 x lO-^i/iioo eV, 3.7 x IQ-^i/iioo eV, and 
5.4 X 10~"^^/iioo eV, respectively. The values of the other 
parameters are the same as those for Fig[T2| except 
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Fig. 13 Spatial power spectrum of the CMB temperature 
anisotropy computed for m^ = 2.0 x 10~^'^/iioo eV, 3.7 x 
10~^^/iioo eV, and 5.4 x lO^^^/iioo eV using our model. The 
values of the remaining parameters are the same as those 
for Fig. 1121 except that fi^^o ~ 0.15 is employed here. The 
WMAP data are also shown for the purpose of comparison. 
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that we adopt m^ = 3.7 x lO^'^^/iioo eV as the most 
appropriate value. 

Figure [T3] shows the CMB power spectra computed 
for f}m,o = 0.182, 0.237, and 0.292. The values of the 
other parameters are the same as those for Fig. [T^l ex- 
cept that nicf, = 3.7 x 10~'^^/iioo eV is employed here. 
The WMAP 3-yr data are also plotted for the purpose 
of comparison. Note that as the value of ^m.o increases, 
the amplitudes of the peaks, particularly the first and 
second peaks, of the spectrum tend to decrease signifi- 
cantly. In other words, the value of rim.o has a strong 
influence on the heights of the peaks, although it hardly 
affects the spatial periodicity and the amplitudes of the 
number count variation of the galaxy N-z relation. We 
select fim,o = 0-237 as the best value based on the com- 
parison of the theoretical curves with the WMAP 3-yr 
data. 

Figure [15] shows the CMB power spectra computed 
for f^b,o^?oo = 0-020, 0.022, 0.024, and 0.026. The val- 
ues of the other parameters are the same as those em- 
ployed for FigHH except that fini.o — 0.237 is adopted 
here. The WMAP 3-yr data are also plotted for com- 
parison. As the value of fib.o increases, the height of 
the first peak of the theoretical spectrum tends to in- 
crease, whereas that of the second peak tends to de- 
crease. Therefore, we select 0-024 for r^b.o^inn as the 
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Fig. 14 Spatial power spectrum of the CMB tempera- 
ture anisotropy computed for rim.o = 0.182, 0.237, and 
0.292 using our model. The values of the remaining pa- 
rameters are the same as those for Fig. 1131 except that 
rriri, = 3.7 x 10~^/iioo eV is employed here. For the purpose 
of comparison, the WMAP 3-yr data are also plotted. 



that n^fi = 0.15 is employed here. For the purpose of 
comparison, the WMAP 3-yr data are also plotted in 
the diagram. We find it rather difficult to select the best 
value for m^ based solely on the CMB power spectrum. 
However, the computed N-z relation depends strongly 
on the value of this parameter, as shown in Fig. El so 
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Fig. 15 Spatial power spectrum of the CMB temperature 
anisotropy computed for Qh,ohwo = 0.020, 0.022, 0.024, 
and 0.026 using our model. The values of the remaining 
parameters are the same as those for Fig. 1141 except that 
i^m.o ~ 0.237 is employed here. For the purpose of compar- 
ison, the WMAP 3-yr data are also plotted. 



In a manner similar to that described above, we 
eventually succeeded in producing a model with a non- 
minimally coupled scalar field that can account not 
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only for the observed spatial periodicity exhibited by 
the N-z relation of the field galaxies but also for the 
observed spatial power spectrum of the CMB tempera- 
ture anisotropy. While this model may not be the true 
best-fit in view of the procedure we have employed to 
deduce the individual values of the model parameters, 
we have established the possibility that the observed 
picket-fence struct ure of the galaxy N-z rel ation origi- 
nally discovered bv iBroadhurst et alJ (|l990r ). which re- 
mains clearly observable in the 2dF data, may be real, 
rather than an observational fluke. A more elaborate 
parameter search to find a better model or models using 
Markov chain Monte Carlo (MCMC) methods (Lewis 
and Bridle 2002) is currently underway. 

3.3 Type la supernovae 

Another important type of observational data that 
should be considered is the Hubble diagram of the la- 
type supernovae (SNI a) based on the Su pernovae Gold 
Dataset compiled by iRiess et al.l (|20Q7l ) , containing a 
total of 182 objects (119 data points fr om the previous 



sample compiled bv lRiess et al.l (|2004l) . 16 data points 



discovered recently using the Hubble Space Telescope 
(HST), and 47 of 73 distant objects included in the first 
year relea se of the Sii pernova Legacy Survey (SNLS) 
sample of lAstier et al.l (2006). We would therefore like 
to test the viability of our best-fit model against this 
diagram. 

Since we are concerned only with the flat geometry 
case, we have the following expression for the luminosity 
distance Z3l in units of the present-day Hubble radius 
c/Ho: 



Dl = {1 + z) I dz 
Jo 



H{z)' 



where H{z) is the Hubble parameter. 
H{z) = {da/dt)/a is 



(21) 



From Eq. lM]) . 



H{z) 



6^H' + \ mnT + (1 - 6C0') 



^')2^fm,c^ 



' exp (- 



+Fo'("m,o(l + zf + r!,.o(l + zf 



1/2 



+ ^A,o) 



(1-6^02), 



(22) 



where a primed quantity designates its derivative with 
respect to the proper time t, viz., (jj' = d(f)/dt. 
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Fig. 16 (a) Hubble diagram or the distance modulus (ordi- 
nate) vs. redshift (abscissa) relation of SNIa's. The discrete 
points with error bars indicate the obse rvational data take n 
from the Supernovae Gold Dataset of I Riess et al.l (|2007J). 
The solid(red) line curve and the broken(green) line curve 
denote the theoretical computations obt ained for ou r scala r 
field model and for the ylCDM model of lRiess et ai] (|2007l '). 
respectively, (b) Computed distance moduli and observa- 
tional data, with the corresponding values given for the 
constant expansion model being subtracted, plotted as func- 
tions of redshift z for better legibility. (Slm.o = 0, fi^.o = 0, 
and JIk.o = —{c/Ho)^K = 1 corresponding to a negative 
value of K. 



Using this luminosity distance Dl, 
distance modulus /i: 



^I: 



M = 5 logi 



Ho 



-Dj 



25 



we obtain the 



(23) 



provided that c/ Hq is given in units of Mpc. 

In Figure I16r a) , we show the theoretical curve of ^ 
obtained for our model (solid (red) line) and that ob- 
tained for the fiat vlCDM model with Q^_o = 0.29 of 
Riess et al.l ( 20071 1 (broken(green) line) as functions of 
redshift z in the range [0.001,2] together with the 182 
data points taken from the Supernova la Gold Dataset 
( Riess et al.l 12007). The values of the model param- 
eters employed for our model are the same as those 
for the best-fit model obtained in the previous section: 
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viz., ^11,0 = 0.237, r20,o = 0.15 (hence, flA.o = 0.613), 
^ = -40, m^ = 3.7 X lO^^i/iioo eV, and q = 12000 
together with /iioo = 0.72 (assumed). Our theoreti- 



the c urvature scalar, as was done bv iMorikawal (1990, 
19911 ). This enables us to account for the picket-fence 



cal cu rve is indistinguishable from that of iRiess et al. 
( 20071 ). In other words , our model is at least as good 



as the ylCDM model of JRiess et~aLl (|2007t ) so far as the 
Hubble diagram of SNla's is concerned. 

In FigurefTHTb). the theoretical distance moduh ^ for 
o ur model(sol i d(red) line curve) and the ylCDM model 
of iRiess et al.l (|2007l ) (broken(green) line curve), as well 



as the observational data, with the corresponding val- 
ues computed for the constant expansion model uni- 
verse (r^nio — 0, fli-,0 = 0, flA,o = 0, and Ok,o = 
— {c/Ho)^K — 1 corresponding to a negative value of 
K) being subtracted, are shown as functions of redshift 
z in the range [0.001,2] for better legibility. Since the 
oscillation of ((> arises only for z < 1 in our model (see 
Fig. [2]), our result differs only slightly from that of the 
ylCDM model for z ~ 1 or greater. The x^ value for 
our model and that for the vlCDM model are 159.4 and 
160.6, respectively. However, the difference between the 
two theoretical curves begins to appear at z less than 
~ 0.02, and the theoretical curve generated by the pro- 
posed model exhibits an oscillatory variation with de- 
creasing z. Hence, it is likely that future observations of 
nearby SNla's may r eveal which type of model is more 
plausible. Note that iLazkoz et a l. (2005) examined 11 
types of cosmological models by comparing the result- 
ing theoretical Hubble diag rams with those o f a num- 
ber of SNla's pubhshed bv iRiess et all (|2004h . to find 
a model employing the equation of state in which the 
variables oscillate with time to be the most adequate. 



4 Conclusion 

We have succeeded in developing a cosmological model 
that can explain, in a unified manner and yet within 
the framework of the flat Robertson- Walker metrics, 
three types of major cosmological observations, viz., 
the picket-fence structure or the spatial periodicity of 
galaxy N-z relation, the spatial power spectrum of the 
CMB temperature anisotropy, and the Hubble diagram 
of recently observed SN la's. The picket-fence st ructure 
originally discovered bv lBroadhurst et al.l (|l990l in the 
N-z relation of the field galaxies observed through their 
penci l-beam survey rema ins distinctly visible in the 2dF 
data (jColless et al.ll200l[ ). as well as in the SDSS data, 
although it is less pronounced. Such a feature is diffi- 
cult to explain with the Robertson- Walker metric cos- 
mology. An ingenious method by which to account for 
the presence of the picket-fence structure is to intro- 
duce a scalar field that is non-minimally coupled with 



structure as an apparent or illusory effect caused by the 
oscillation of the expansion rate with respect to time. 
We have therefore investigated if it would be possible 
to produce a model, accompanied by Morikawa's scalar 
field, that is capable of accounting simultaneously for 
the picket-fence structure of the N-z relation obtained 
from the 2dF data and the spatial power spectrum of 
the CMB temperature anisotropy based on the WMAP 
3-yr data. However, the best-fit model deduced from 
the 2dF data and that derived from the WMAP data 
have turned out to be mutually exclusive so long as we 
adhere to the strict Morikawa scalar field. 

The most straightforward method by which to avoid 
this difficulty was to modify the form of the potential 
for the scalar field. Based on a number of experiments, 
we have ventured to introduce a new potential func- 
tion that is proportional to (jP exp(— qc^^), with g being 
constant instead of <jy^ employed bv IMorikawal (1990, 
19911 ) . The effect of q is to retard the time at which 



starts growing, moving it closer to the present epoch. 
Furthermore, we need to maintain the scalar field at 
an extremely small magnitude throughout the epochs 
from the birth of the CMB to the onset of the growth 
of (f). In addition, we have found that by solving the 
evolution equations in the forward direction, starting 
with such a stationary state at a = 10~^ and a suit- 
ably chosen value for the power index q, we can pro- 
duce a cosmological model in which (f> possesses the 
desired behavior without requiring any particular fine- 
tuning of the initial condition. For q = 12, 000, we are 
now able to reproduce the observed power spectrum of 
the CMB temperature anisotropy as well as the picket- 
fence structure exhibited by the A'^-z relation by using 
a model that is characterized by the following values of 
the model parameters: flm.o — 0.237, r2b,o'iioo — 0.024, 



n 



</),0 



0.15(hence fl 



0.613), C = -40, TO0 



A,0 

f>o = 2.5206 X lO^'^ (corresponding 
init == 6.417 X 10-31), and n^ = 1 



3.7 X lO-31/iioo eV 
to the initial value 
(Harrison-Zel'dovich spectrum) together with the as- 
sumption that hiQo = 0.72 and rcMB,o = 2.726 K. 

Figure fTW a) shows the time evolution of the scale 
factor of the proposed m odel in compa ri son w ith that 
of the ylCDM model of ISpergel et al! (|2007t ). Note 
that the rate of expansion of the proposed model is 
somewhat lower a t Hot above ~ 3, c ompared with the 



ylCDM model of lSpergel et al.l ([20071) due to the smaller 
value adopted for the cosmological constant A in the 
proposed model. Figure [TTfb) is an enlargement of the 
area of Fig. fTTT a) in the range of [0.4, 0.6] for Hgt, which 
shows a wavy behavior of a for the proposed model be- 
cause the scalar field is in the phase of a damping os- 
cillation here. The present age of our model universe 
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turns out to be 13.2 Gyrs, whic h is somewhat lowe r 
than that of the ylCDM model of ISpergel et alj (|2007t) . 
Neverthel ess, the lower l imit f or the cosmic age of 13 
Gy set by iKomiva et alJ (|2006l ) is cleared. 




Fig. 17 (a) Comparison of the time evolution of the 
scale factor of our model with that of the /ICDM model of 
ISpergol ot al. (2007). The abscissa shows Hot, the proper 
time t normalized to the Hubble time, (b) Enlargement of 
the part of Fig. [l^a) in the range of [0.4, 0.6] for Hot. 

The proposed model has also been tested against the 
observed Hubble diagram of the S NIa's created from 

()2007l) . 



the Supernovae la Gold Dataset of iRiess et al 



We have found that the proposed i nodel is at l east a s 
satisfactory as the /ICDM model of iRiess et al.l (|2007l ). 
although these two models show a rather noticeable dif- 
ference for z less than approximately 0.02. 

We would now like to turn our attention to the ef- 
fect of the possible change of the gravitational con- 
stant G ( s ee Eq . iirSI) ) on the luminosities of galaxies. 
Hill et al.l ( 1990h . for instance, suggested that the spa- 



tia l periodicity of the gala xy N-z relation discovered 
by iBroadhurst et al.l (|l990l ) is due primarily to a peri- 
odic change in the luminosities of galaxies, which might 
be caused by the variation of the gravitational con- 
stant. We have therefore carried out sample compu- 
tations of the N-z relation, taking into account the 
gravitational correction (G-correction) Eq{z) in obtain- 
ing the apparent magnitudes of galaxies employing the 
proposed best-fit m odel deduced in th e prev ious sec- 
tions. According to [Hansen fc Kawaleij (|l994l ). the lu- 
minosity L of a galaxy is proportional to G^"^*, so 
that we employ the fol lowing form for the G-correction 



( Kawabataet al.ll2002f ): 



Eciz) = -2.5log 



GcS,o 



7.5/3 




0.07 



0.08 



0.09 



0.12 



Fig. 18 Effect of G-correction on the N-z relation of galax- 
ies. The abscissa is the redshift z, and the ordinate shows 
the fractional counts A'^, as in the previous figures. The 
solid(red) line curve shows the number counts of galaxies 
taking into consideration the G-correction arising from the 
temporal variation of the gravitational constant G, while 
the broken(green) line curve indicates the result without G- 
correction. The hatched area shows the observational data 
of the 2dF survey. 



where Gcff,o is the effective value of the gravitational 
constant at the present era, while /? is a positive fudge 
factor (< 1) to account for the uncertainty involved 
in the power raised to G. For simplicity, however, we 
adopt f3 — I for our computations. 

The solid(red) line curve and the broken(green) line 
curve in Fig. [THl denote the resulting N-z relations for 
the range 0.07 < z < 0.12 obtained with and with- 
out the G-correction, respectively, whereas the hatched 
area shows the 2dF data. At z ~ 0.1, for instance, 
the value of G^s/G is no less than 0.99, as shown in 
Fig. ^d) (see also Eq.([TH])). Equation ^ predicts 
that Eq{z) shall not exceed 0.08 in magnitude, which 
in turn indicates that the change in the number counts 
N caused by the temporal variation of the gravitational 
constant G should be less than 0.61%. Hence, we may 
safely conclude that the primary cause for the picket- 
fence structure of the spatial periodicity of the galaxy 
distribution is the oscillatory expansion rate of the uni- 
verse arising from the oscillation of the non-minimally 
coupled scalar field, and that the astrophysical effect 
originating from the variation of G must be at most 
secondary. 

Finally, the value of the Brans-Dicke parameter wbd 
for our model is found to be 



(24) 



wbd 



1 - 6$0§ 
2Ae4>l 



40, 986 



(25) 



16 



which is not inconsistent with the lower hmit ^ 40, 000, 
referred to as the Cassini bound, which w as set based on 



recent observations of our solar system (jBertotti et al. 
[2003; .Berti et al.:.2005) . 
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A The Perturbation Equations 

In this Appendix, we present the formulation fo r the e v olution of scalar perturb a tions of our model. We closely 
follow the me thod employed by Baccigalupi et alj (|2000l ): IPerrotta fc Baccigalupil (|2002l ) based on the formalism 
developed bv lMa fc Bertschingeij (jl995l ) to describe the evolution of perturbations in the synchronous gauge. In 
this type of gauge, the perturbed metric takes the form 



where hij represents the metric perturbations, and its Fourier transform can be written as 

1 



i^ij \^-) ^c) 



d^ke'^'' 



kikjft,(k, tc 



kikj - -5ij ] 6ri{k,tc 



(Al) 



(A2) 



where h denotes the trace of hij, and ry represents its traceless component. 

As for the perturbation c omputations, we employ the formulation reported by iBaccigalupi et al.l ( 20001 ); 
Perrotta fc Baccigalupil (|2002l ). By setting 



F = 



8nG 



^r 



(A3) 



^=2 



1 / m^c 



tp exp 



47rG 



■V' 



(A4) 



in their formulae, we obtain our perturbation equations. 

The perturbed densit y (5/3, pressure SP, velocity dive rgence 9, and shear a obtained based on the formalism of 
Baccigalupi et al.l (|2000l ): IPerrotta fc Baccigalupil (J2002l ) are as follows: 



Spc^ = 



1 



1 



M' 



1 - SttG^VVc^ 

^ 2 - inG 



a^ a^ a'^ a'^ 



AttG 



+ (!!^) ^5^(1 _ ,|^^2) ,^p(_g|lj^2) ^ 2^^^^^ ^ ^^^-^ 



(A5) 



6P 



1 



SPrr. 



1 - SttG^V'V"* 
-2^-^S^pip - 2^^i;S^ - A^^^Sip 

tO ^O /-,Z 






L^S'i/j - Ai-^^6^P + 2^'^ijSiP - 2^\s^Pi^ - 2^ij^Si^ 



a 
4fc2 

3 a' 



^iip5%l)' 



w 

2 h 

3 a 



\ h 



AttG 



AttG 



^5ip{\ - q—^ip )exp(-q— ^V ) 



3c4 



3c4 



^itp^ 



(A6) 



(pc^ + p)e 



1 - SttG^VVc^ 



(Pr, 



P, 



va,r j^m.r 



k'' ( - - ■ - a - 

o^ a 



(A7) 



[pc^ + P)a 



1 - STrGCV-VC 



2 fc^ I - 3 - - h 

{pm,rC' + Prn^r)crrn.r + ^—{ -2^^^' + T^ (-2^l/iV^) (?] + -) 



3a2 



6' 



(A8) 



where m and r indicate matter and radiation, respectively. Sip is the perturbed part of the scalar field, and ip 
is the unperturbed part of the scalar field. 



^p = jp -\- Sip. 



(A9) 
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The perturbed Klein-Gordon equation is 

S^ + 2-S^p + k^S^p + 6^-S^p + a^^tpSR + a'^(d^d^V)5i: + -,/; = 0. (AlO) 

a a 1- -r 2 

Using these equations, the perturbed Einstein equations are 

^h-2k^i^^^aHpc^, (All) 



2k^^^aH-pc^ + P)e, (A12) 

i(2fc2^-/i-2-/i) = ^a2jP, (A13) 



-\(h + 6f/ + 2-(/i + 6?7) - 2fc2r,) = ^flS^pc^ + P)a. (A14) 

Using these perturbed Einstein equations, we calculate the perturbations on the theoretical CMB anisotropy 
spectrum. 
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